Kitaev-Heisenberg Model on a Honeycomb Lattice: Possible Exotic Phases in Iridium 
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We derive and study a spin one-half Hamiltonian on a honeycomb lattice describing the exchange 
interactions between Ir 4+ ions in a family of layered iridates ^IrOa (A=~Li, Na). Depending on the 
microscopic parameters, the Hamiltonian interpolates between the Heisenberg and exactly solvable 
Kitaev models. Exact diagonalization and a complementary spin-wave analysis reveal the presence 
of an extended spin-liquid phase near the Kitaev limit and a conventional Neel state close to the 
Heisenberg limit. The two phases are separated by an unusual stripy antiferromagnetic state, which 
is the exact ground state of the model at the midpoint between two limits. 

PACS numbers: 75.10.Jm, 75.25.Dk, 75.30.Et 



Magnetic systems exhibit, most commonly, long-range 
classical order at sufficiently low temperatures. An ex- 
ception are frustrated magnets, in which the topology of 
the underlying lattice and/or competing interactions lead 
to an extensively degenerate manifold of classical states. 
In such systems, exotic quantum phases of Mott insula- 
tors (spin liquids, valence bond solids, etc.) can emerge 
as the true ground states (for reviews see Refs. 0,0). In 
quantum spin liquids, strong zero-point fluctuations of 
correlated spins prevent them to "freeze" into magnetic 
or statically dimerized patterns, and conventional phase 
transitions that break time-reversal and lattice symme- 
tries are avoided. Spin liquids have attracted particular 
attention since Anderson proposed their possible connec- 
tion to superconductivity of cuprates 

Recently, spin- liquid states of matter have been exem- 
plified, on a quantitative level, by an exactly solvable 
model by Kitaev Q ■ His model deals with spins one-half 
that live on a honeycomb lattice. The nearest-neighbor 
(NN) spins interact in a simple Ising-like fashion but, be- 
cause different bonds use different spin components [see 
Fig-da)], the model is highly frustrated. Its ground state 
is spin-disordered and supports the emergent gapless ex- 
citations represented by Majorana fermions (JJ. Spin- 
spin correlations are, however, short-ranged and confined 
to NN pairs 0, @] . This may suggest the robustness of 
the disordered state to spin perturbations. Indeed, Tsve- 
lik has shown [?[ that there is a window of stability for 
the spin-liquid state in the Kitaev model perturbed by 
isotropic Heisenberg exchange. 

Finding a physical realization of this remarkable model 
is a great challenge, also because of its special properties 
attractive for quantum computation [4]. As the key el- 
ement of the model is a bond-selective spin anisotropy, 
one possible idea is to explore Mott insulators of late 
transition metal ions with orbital degeneracy, in which 
the bond directional nature of electron orbitals can be 
translated into a desired anisotropy of magnetic interac- 
tions through strong spin-orbit coupling. 



In this Letter, we examine the iridium oxides A2ITO3 
from this perspective. In these compounds, the Ir 4+ ions 
have an effective spin one-half moment and form weakly 
coupled honeycomb-lattice planes. Our analysis of the 
underlying exchange mechanisms shows that the spin 
Hamiltonian comprises two terms, ferromagnetic (FM) 
and antiferromagnetic (AF), in the form of Kitaev and 
Heisenberg models, respectively. The model has an in- 
teresting phase behavior and hosts, in addition to the 
spin-liquid state, an unusual AF order that is also an 
exact solution at a certain point in phase space. 

Experimental studies of iridium compounds are rather 
scarce, and the nature of their insulating behavior is not 
yet fully understood. In fact, Na2lr03 was suggested as 
an interesting candidate for a topological band insulator 
Given that high temperature magnetic susceptibil- 
ities of Na2lr03 and Li2lrC>3 obey the Curie- Weiss law 
wit h an ef fective moment corresponding S = 1/2 per Ir 



13" HI) we start here with the Mott insulator picture. 
The Hamiltonian.- We recall that the Ir 4+ ion in the 
octahedral field has a single hole in the threefold degener- 
ate t2g level hosting an orbital angular momentum 1 = 1. 
Strong spin-orbit coupling lifts this degeneracy, and the 
resulting ground state is a Kramers doublet with total 
angular momentum one-half [3], referred to as "spin" 
hereafter. In fact, it is predominantly of orbital origin, 
and this is what makes the magnetic interactions highly 
anisotropic due to the spin-orbit entanglement of mag- 
netic and real spaces. In ^4 2 Ir03 compounds, the Ir06 
octahedra share the edges, and Ir ions can communicate 
through two 90° Ir-O-Ir exchange paths Q or via direct 
overlap of their orbitals. Collecting the possible exchange 
processes (discussed below) and projecting them onto the 
lowest Kramers doublet with S = 1/2, we obtain the fol- 
lowing spin Hamiltonian on a given NN ij bond: 
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(1) 



Here, spin quantization axes are taken along the cubic 
axes of IrOg octahedra. In a honeycomb lattice formed 
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by Ir ions, there are three distinct types of NN bonds 
referred to as r y(= x, y, z) bonds because they host the 
Ising-like J\ coupling between the 7 components of spins 
[see Fig. QJa)]. The first part of Eq. ([T]) is thus nothing 
but the FM Kitaev model, and the J2 term is a conven- 
tional AF Heisenberg model. The exchange constants J\ 
and J2 are derived from a multiorbital Hubbard Hamilto- 
nian consisting of the local interactions and the hopping 
term. The latter describes t p d-K hopping between Ir 5d 
and O 2p orbitals via the charge-transfer gap A p d, and 
a direct dd overlap t' between NN Ir t% g orbitals [lij ]. 
We find J\ — [r\\ + 2772) and J2 = (772 + Hereafter, 
we use 4t 2 /9Ud as our energy unit, where t = t pd ^/A p d, 
and Ud stands for the Coulomb repulsion on the same 
d orbitals. There are three physically distinct virtual 
processes that determine the set of r\ parameters and 
thus the ratio Jij 3\. The 7/1 = v ^" Jh u^j h term ap- 
pears due to the multiplet structure of the excited lev- 
els induced by Hund's coupling Jh 0- The processes 
when two holes meet at the same oxygen site (and ex- 
perience Up repulsion) and when they are cyclically ex- 
changed around a 1x2 O2 square plaquette bring together 
a 772 = a — U , f tt in -^ jL - contribution. Further, a direct dd- 
hopping t' between NN Ir ti g orbitals contributes to the 
Heisenberg term with exchange coupling 773 = {t 1 /t) 2 . It 
is difficult to estimate the values of all the parameters 
involved; however, we expect 771 to be the largest, of the 
order of 1, and 772,3 < 1. 

We parametrize the exchange couplings as J\ = 2a and 
J2 = 1 — a and study the properties of Kitaev- Heisenberg 
model ([1]) in the whole parameter space < a < 1. 

Phase diagram.- At a — 0, we are left with the Heisen- 
berg model exhibiting the Neel order with a staggered 
moment reduced to {S z } ~ 0.24 [l6|. The opposite limit, 
a = 1, corresponds to the exactly solvable Kitaev model 
with a short-range spin- liquid state Q , where spin corre- 
lation functions are identically zero beyond the NN dis- 
tance and, on a given NN bond, only the components of 
spins matching the bond type are correlated [5j. 

Interestingly, the model is exactly solvable at a = h, 
too. At this point Eq. ([1} reads, e.g., on a z-type bond, 



(a) 



as U\f = \{S'fS* + SfS] - SfS?). This anisotropic 
Hamiltonian can be mapped to that of a simple Heisen- 
berg model on all bonds simultaneously [I?} • Specifically, 
we divide the honeycomb lattice into four sublattices [see 
Fig. Hfb)] and introduce the rotated operators S: While 
S = S in one of the sublattices, S on the remaining 
three sublattices differs from the original S by the sign 
of two appropriate components, depending on the sub- 
lattice they belong to. In the new basis, Eq. ([T]) takes 
the form 



W|7> = -2(2a 



1) S?S] -(1- a) Si-Sj 



(2) 



At a 



i, the first term vanishes and we obtain the 
isotropic, both in spin and real spaces, Heisenberg model 
' - \ Si-Sj with FM coupling. Thus, at a = |, 
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FIG. 1: (a) Three types of bonds in the honeycomb lat- 
tice and Kitaev part of the interaction, (b) The supercell 
of the four-sublattice system enabling the transformation of 
the model ([T]) into the Hamiltonian of a simple ferromagnet 
at a = i. This supercell with periodic boundary conditions 
applied was used as a cluster for the exact diagonalization. 
(c) Schematic phase diagram: With increasing a, the ground 
state changes from the Neel AF order to the stripy AF state 
(being a fluctuation- free exact solution at a = |) and to the 
Kitaev spin liquid. See the text for the critical values of a. 



i.e., at Ji = 2J2, the exact ground state of model ([T]) 
is a fully polarized FM state in the rotated basis. Now 
consider the FM array of spins with, e.g., (S z ) = 1/2, 
and map it back to the original spin basis. The resulting 
order corresponds to a stripy AF pattern of the original 
magnetic moments depicted in Fig. QJc). Note that such 
a stripy order, despite being of AF type, is fluctuation- 
free at a = h and would thus show a fully saturated AF 
order parameter. 

The above discussion suggests three possible ground 
state phases of the model |T|) as shown in Fig. HJc): (i) 
Neel order near a — 0, (ii) stripy AF order around a — ^, 
and (iii) a spin-liquid phase close to a = 1. 

We first consider the ordered phases. Except special 
cases of a — and a — i just discussed, the Hamil- 
tonian ([]]) does not have any spin-rotational symmetry. 
However, a spurious SU(2) continuous symmetry and as- 
sociated pseudo-Goldstone mode appear in a linear spin- 
wave (SW) descri ptio n. As in the case of a similar model 
on a cubic lattice |18| , we find that quantum fluctuations 
restore the underlying discrete (hexagonal) symmetry of 
the model, selecting thereby the direction of ordered mo- 
ments along one of the cubic axes (of IrOg octahedra), 
and also open a gap in SW spectra. Considering the 
quantum energy cost for rotating the order parameter by 
a small angle away from a cubic axis, we find a quantum 
SWgapA~f(a-i) 2 fora~i. 
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The classical phase boundary between Neel and stripy 
AF orderings is at a = i, where linear SW spectra of 
both states develop zero-energy lines [ljij], reflecting the 
infinite degeneracy of classical states. At a = |, Eq. (fTJ) 

reads, e.g., on z-type bonds, as H$ = §(Sf Sf + S?SJ); 
i.e., only two spin components are coupled on a given 
bond. Considering Neel or stripy AF with ordered spins 
parallel to the z axis, one finds that flipping all the spins 
along a zig zag chain, formed by x-and y-type bonds, does 
not change classical energy. This degeneracy is again ac- 
cidental (an artifact of classical treatment) and can thus 
be lifted by quantum fluctuations. They favor the Neel 
state and shift the classical phase boundary to a larger 
value a ~ 0.4. This estimate is obtained by comparing 
the energies of the Neel [ei ~ —^(3 — 5a)] and the stripy 
[e2 — — f (5a — 3 + —)] states including quantum correc- 
tions via second-order perturbation theory and matches 
well the numerical result found below. 

Now, we discuss the phase behavior at | < a < 1, 
i.e., in between two exact solutions (stripy AF at a = ^ 
and a Kitaev spin- liquid at a = 1). In terms of rotated 
spins, all the couplings are of FM nature in this region 
[see Eq. (J5J]. Thus, the FM order (read stripy AF of the 
original spins) is the only possible magnetic phase here 
to compete with the spin-liquid state. Since the latter is 
stable against a weak Heisenberg-type perturbation Q, 
a critical value of a for the spin order/disorder transi- 
tion must be located at some point less than 1. We give 
its naive estimate based on the energetics of these two 
phases. The energy of the stripy AF state is given above. 
The upper boundary for the energy of spin- liquid state is 
given by the expectation value of Eq. ^ using the exact 
result (S^Sj) — 0.13 at a = 1 [5J. As a result, we find 
the transition from stripy AF order to a spin liquid at 
a ~ 0.86 (close to the numerical result below). 

Single-magnon excitations fail to detect this transition 
(since, as said above, there is not any other competing 
magnetic state). As a increases, the lower branch of the 
linear SW spectrum just gradually softens, to become 
completely flat in the limit of a = 1 where the classical 
ground state is extensively degenerate [13] • We therefore 
suspect that the instability responsible for the collapse 
of magnetic order resides in the two-magnon sector [21| . 
Leaving this subtle issue for a future work, we now turn 
to our numerical results, which describe the evolution of 
spin correlations across the entire phase diagram. 

Numerical study - We use the Lanczos exact diagonal- 
ization method to study a 24-site cluster [see Fig. QJb)] 
with periodic boundary conditions. The cluster is com- 
patible with the above discussed four-sublattice transfor- 
mation of Eq. (fTJ into Eq. ([2]). This provides an exact 
reference point a = 1/2, which is useful for the interpre- 
tation of numerical data shown in Figs. ^ and ([3]) in 
terms of the original as well as transformed spins. 

Figure [2] clearly locates the two phase transitions. In 
particular, a pronounced maximum in the second deriva- 
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FIG. 2: (a) Squared total spin of the 24-site cluster, normal- 
ized to its value in the fully polarized FM state, as a function 
of q. The solid (dot-dashed) line corresponds to the rotated 
(original) spin basis, (b) The NN spin correlations: The solid 
(dot-dashed) line corresponds to a scalar product of the ro- 
tated (original) spins. The component of the correlation func- 
tion matching the bond direction is indicated by a dotted line. 
This quantity is the same in both bases. The inset compares 
NN spin correlations (solid line) above a — 0.5 with longer 
range spin correlations up to third-nearest neighbors (dotted 
lines), (c) Negatively taken second derivative of the ground 
state energy with respect to a. Its maxima indicate the phase 
transitions at a ~ 0.4 and 0.8. 



tive of the ground state energy [Fig.^c)] indicates a first- 
order transition from Neel to stripy AF phase at a ~ 0.4. 
The much weaker (note the log scale) and wider second 
peak at a ~ 0.8 suggests a second- (or a weakly first-) 
order transition from stripy AF to a spin-liquid state. 

Figure[5]Ja) shows the squared total spins sl ot and S^ ot 
normalized to S(S+1) with S = N/2 that can be reached 
in the FM state. Although these are not conserved quan- 
tities in the model, they characterize the phase map quite 

-2 

well. In particular, a long tail of S tot above a = 0.8 in- 
dicates a "leakage" of stripy AF correlations into a spin- 
liquid phase. This is also evidenced by the behavior of 
longer range, beyond NN, spin correlations that are still 
visible in a spin-liquid regime, except close to the Ki- 
taev limit where they vanish completely [see the inset in 
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FIG. 3: Magnetic moment 2(S t z ot )/iV induced by field B z 
(Zeeman coupled to the rotated spins). The circles at a > 0.8 
show the inverse spin susceptibility to this field. 



in the Kitaev model as well 23 1. 

In conclusion, we have examined the interactions and 
possible magnetic states in iridates A 2 Ir03. The ob- 
tained Kitaev-Heisenberg model shows rich behavior in- 
cluding a spin liquid and unusual stripy AF phases. We 
hope that these results will motivate experimental studies 
of layered iridates and similar compounds of late transi- 
tion metal ions, where the physics of the Kitaev model 
might be within reach. 

We thank B. Keimer, A. Schnyder, S. Trebst, and M. 
Zhitomirsky for discussions. We are grateful to H. Tak- 
agi and A.M. Tsvelik for valuable discussions and for 
communicating the unpublished results. Support from 
MSM0021622410 (J.C.) and GNSF/ST09-447 (G.J.) is 
acknowledged. 



Fig. I^b)]. 

Figure HJb) highlights how the NN spin correlations 
evolve as their interactions change from one type to 
another. In the Neel state, where the model is more 
Heisenberg-like for the original spins, we reproduce 
(SrSj) ~ -0.37 0. At the "hidden" FM Heisenberg 
point a = 1/2, one finds (Si ■ Sj) = j, equally con- 
tributed by all three components of the rotated spin 
S. Things change dramatically in the spin-liquid phase: 
Here, a particular component of spin correlations (SjiSj), 
dictated by the Kitaev model, dominates. Its value of 
0.132 that we find at a = 1 agrees well with the exact 
result 0.131 for an infinite lattice [f|. 

Finally, we discuss the response to a weak magnetic 
field B z which, in terms of original spins, linearly couples 
to the stripy AF order parameter. Figure [3] shows that 
even a very weak field induces a nearly saturated moment 
in the entire region of the stripy AF phase. As the sys- 
tem switches to the Neel phase, a response to the "stripy 
field" B z drops abruptly to zero, as expected. The in- 
duced moment sharply reduces near a — 0.8, too, but 
remains finite in a spin-liquid phase. Here the magnetiza- 
tion curve shows a linear dependence on B z , and we may 
extract from its slope the susceptibility \ = (^tot)/^^- 
Shown in Fig. [3] is the inverse value of x as a function of 
a. This quantity scales with the energy gap between the 
ground state and the excited states that are accessible 
by the magnetic field. According to Kitaev's solution 
these states must belong to the flux sectors located at 
energies of the order of 1. The observed x 1 oc (cv — 0.8) 
behavior shows that this characteristic (spin) gap grad- 
ually softens towards the a ~ 0.8 critical point, as the 
spin correlations beyond the NN distances start to grow 
[see Fig. (Kb), inset]. 

Experimental data are rather insufficient to 

conclusively locate the position of A2L-O3 comjpounds in 
our phase diagram. Also, Na/Ir site disorder [13j has to 
be kept in mind: Often, nonmagnetic impurities induce 
local moments [22], and this has been shown to happen 
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